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Abstract 

The Homotopy Perturbation Method (HPM) is used to solve the Burgers-Huxley non-linear diffe- 
rential equations. Three case study problems of Burgers-Huxley are solved using the HPM and the 
exact solutions are obtained. The rapid convergence towards the exact solutions of HPM is nu- 
merically shown. Results show that the HPM is efficient method with acceptable accuracy to solve 
the Burgers-Huxley equation. Also, the results prove that the method is an efficient and powerful 
algorithm to construct the exact solution of non-linear differential equations. 
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1. Introduction 

Most of the nonlinear differential equations do not have an analytical solution. Recently, semi-analytical so- 
lutions of real-life mathematical modeling are considered as a key tool to solve nonlinear differential equa- 
tions. 

The idea of the Homotopy Perturbation Method (HPM) which is a semi-analytical method was first pioneered 
by He [1], Later, the method is applied by He [2] to solve the non-linear non-homogeneous partial differential 
equations. Nourazar et al. [3] used the homotopy perturbation method to find exact solution of Newell-Whitehead- 
Segel equation. Krisnangkura et al. [4] obtained exact traveling wave solutions of the generalized Burgers- 
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Huxley equation by using hyperbolic tangent method. New exact solutions of the generalized Burgers-Huxley 
equation are also obtained by Gao and Zhao [5], Hashim et al. [6] have solved the generalized Burgers-Huxley 
equation by using adomian decomposition method. Wang et al. [7] studied the solitary wave solutions of the 
generalized Burgers-Huxley equation. Also, Darvishi et al. [8 1 have used spectral collocation method and Dar- 
vishi’s preconditionings to solve the generalized Burgers-Huxley equation. 

The generalized Burgers-Huxley equation is investigated by Satsuma [9] as: 

du d 1 u „ du _ /, „\t „ \ 

— = — — - au — + J3u(l-u ) [u -/), 0<x<l, f>0 (1.1) 

dt dx~ dx v ' v ' 

where a,j3> 0 are real constants and n is a positive integer and y e [0, l] . Equation (1.1) models the inte- 
raction between reaction mechanisms, convection effects and diffusion transports [9], 

When a = 0 and n = 1 , Equation (1.1) is reduced to the Huxley equation which describes nerve pulse pro- 
pagationin nerve fibers and wall motion in liquid crystals [10], When [:’> = 0 and n = 1 , Equation (1.1) is re- 
duced to the Burgers equation describing the far field of wave propagation in nonlinear dissipative systems [11]. 
When n = 1, a ^ 0 and fi ^ 0 , Equation (1.1) is turned into the Burgers-Huxley equation. The Burgers-Hux- 
ley equation is used to model the interaction between reaction mechanisms, convection effects and diffusion 
transport, nerve pulse propagation in nerve fibers as well as wall motion in liquid crystals. 

In the present research work, the Homotopy Perturbation Method (HPM) is applied to obtain the closed form 
solution of the non-linear Burgers-Huxley equation. Three case study problems of non-linear Burgers-Huxley 
equations are solved by using the HPM. The trend of the rapid convergence towards the exact solution is shown 
when compared to the exact solution. 

The idea of homotopy perturbation method is presented in Section 2. Application of the homotopy perturba- 
tion method to the exact solution of Burgers-Huxley equation is presented in Section 3. 

2. The Idea of Homotopy Perturbation Method 

The Homotopy Perturbation Method (HPM) is originally initiated by He [1], This is a combination of the clas- 
sical perturbation technique and homotopy technique. The basic idea of the HPM for solving nonlinear differen- 
tial equations is as follow; consider the following differential equation: 

E(u) = 0, (2.1) 

where E is any differential operator. We construct a homotopy as follow: 

H (u, p) = (l- p)F(u)+ p(e(u)-F(u)} . (2.2) 

where F ( u ) is a functional operator with the known solution v 0 . It is clear that when p is equal to zero 
then //(m,0) = F(m) = 0, and when p is equal to 1, then H (u, \)= E (it ) = 0 . Tt is worth noting that as the 
embedding parameter p increases monotonically from zero to unity the zero order solution v 0 continuously 
deforms into the original problem E(u ) = 0 . The embedding parameter, p e [0, l] , is considered as an expand- 
ing parameter [2], In the homotopy perturbation method the embedding parameter p is used to get series ex- 
pansion for solution as: 

u = YZo P v i = v o + P v i + P\ + p\ + ■ • • ( 2 -3) 

when p — > 1 , then Equation (2.2) becomes the approximate solution to Equation (2.1) as: 

u = v 0 + Vj + v 2 + v 3 H — (2.4) 

The series Equation (2.4) is a convergent series and the rate of convergence depends on the nature of Equation 
(2.1) [1] [2], It is also assumed that Equation (2.2) has a unique solution and by comparing the like powers of 
p the solution of various orders is obtained. These solutions are obtained using the Maple package. 

3. The Burgers-Huxley Equation 

To illustrate the capability and reliability of the method, three cases of nonlinear diffusion equations are pre- 
sented. 
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Case I: in this case we will examine the Burgers-Huxley equation for a = 0, n = 1, y = 1, /? = 1 , so, the equ- 
ation is written as: 

~|2 



du d u , w 



Subject to initial condition: 



<(x,0) = - 



44 

X 

* 4 



44 44 

Jf JC 

4 i « 4 



e ** + e 

We construct a homotopy for Equation (3.1) in the following form: 

du du 0 



H(u,p) = (l-p) 



dt dt 



+ P 



du d'u , w >. 



The solution of Equation (3.1) can be written as a power series in p as: 

u = v Q + pv l + p 2 u 2 H — 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



Substituting Equation (3.4) and Equation (3.2) into Equation (3.3) and equating the terms with identical powers 
of p: 



j0 . dv 0 _ du 0 
dt dt 



>(x,0) = - 



44 

x 

* 4 



44 42 

x x 

4 , „ 4 



1 3 Duq 3 Uq / \ / 1 \ 

p : if + ir“ +u 4 

2 du-, d 2 v, , „ 2 

P ■^- = ^T +Av o u i- u i~ 3u i v o’ 
dt dx 

•5 3^U 0 9 ? 9 / \ 

/7 - -f- + 4y 0 y 2 -y 2 -3y 2 y 0 + 2u x -3u 0 ^ , t> 3 (*,0) = 0. 

dt ox 



+e 

v { (x,0) = 0, 
v 2 (x,0) = 0, 



(3.5) 



Using the Maple package to solve recursive sequences. Equation (3.5), we obtain the followings: 

1 



y 0 (x,t) = - 



V2 

4 A l 

> U i( x ’ t ) = -T 



42 _44 

! 4 +e 4 
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(3.6) 
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By setting p = 1 in Equation (3.4) the solution of Equation (3.1) can be obtained as v = v (t + u, + u, + u 3 + ■ 
Therefore the solution of Equation (3.1) is written as: 



44 



4 42 
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4 



f(x,t) = - 



42 > 
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r .r \ 2 r \ 2 
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(3.7) 
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The Taylor series expansion for 
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is written as: 



77 j_ 

1 1 e 4 ' 4 



a/7 t 
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4 4 
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(3.8) 



Combining Equation (3.8) with Equation (3.7), we get as follow: 

77 



77 . 

JC+- 
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>(x,t) = h 



= — + — tanh 



1 f r 
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(3.9) 



2 2 -A x+ L 2 2 

e 4 4 +e 4 4 

This is the exact solution of the problem. Equation (3.1). Table 1 shows the trend of rapid convergence of the 
results of .S', (xj'j = v 0 (xj'j to S 6 (xj'j = v. (xj'j using the HPM. The rapid convergence of the solution 

toward the exact solution, the maximum relative error of less than 0.0000058% is achieved as shown in Table 1. 
Case II: In Equation (1.1) for a = -1, n = 1, y = 1, p = 1 , the Burgers-Huxley equation is written as: 

du d 2 u du t, w ,> 

— = — - + u — + u(l-u)(u-l) (3.10) 

dt dx 2 dx V A ’ 



Table 1 The percentage of relative errors of the results of (x,t) = v 0 (x,f) to S 6 (.x:,f) = ^ 5 _ of the HPM solu- 

tion of Equation (3.1). 

Percentage of relative error (%RE) 



x — \ x=2 x = 3 



S,(r,r) 


0.01693168743 


0.01002710463 


0.005488150424 


S,(x,t) 


0.000002337346256 


2.025 644856e-7 


9.527700575e-7 


S 5 (x,t) 


2.153484215e-10 


3.464089104e-10 


3.155246333e-10 


S 6 {x,t) 


9.744801271e-12 


3.920421861e-ll 


9.937961084e-ll 


S,{x,t) 


0.05344388963 


0.03164997413 


0.01732302882 


5,(.r,f) 


0.00006806597676 


0.000003967422423 


0.00002580410708 


S 5 (x,t) 


6.184344787e-8 


9.4673 17545e-8 


5.516785201e-8 


S 6 (x,t) 


1 .0087930 18e-8 


1.1 6642282 le-9 


2.026398250e-9 


S,(x,t) 


0.07311570399 


0.04329980772 


0.02369935005 


S,(x,t) 


0.0001675410515 


0.000007498069748 


0.00006123279256 


S 5 {x,t) 


2.799248652e-7 


4.009848425 e-7 


2.364525 175e-7 


S 6 {x,t) 


5.775483086e-8 


6.758498 122e-9 


1.1111 2845 8e— 8 




